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Abstract 

We study the disentanglement of two spin qubits which interact with a general XY spin-chain 
environment. The dynamical process of the disentanglement is numerically and analytically in- 
vestigated in the vicinity of quantum phase transition (QPT) of the spin chain in both weak and 
strong coupling cases. We find that the disentanglement of the two qubits is in general enhanced 
greatly when the environmental spin chain is exposed to QPT. We give a detailed analysis to fa- 
cilitate the understanding of the QPT-enhanced decaying behavior of the disentanglement factor. 
Furthermore, the scaling behavior in the disentanglement dynamics is also revealed and analyzed. 

PACS numbers: 03.65.Vf, 75.10.Pq, 05.30.Pr, 42.50.Vk 
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The coupling between an entangled quantum system and its environment leads to dis- 
entanglement of the system, the process through which quantum information is degraded. 
Disentanglement is a crucial issue that is of fundamental interest due to the fact that the 
distributed nonlocal coherence among multi-particles by the entanglement really matters in 
many important applications of quantum information |2j. Consequently, the fragility of 
nonlocal entanglement is recognized as a main obstacle to realizing quantum computing and 
quantum information processing (QIP) 3|, |4|. Apart from the important link to QIP realiza- 
tions, a deeper understanding of disentanglement is also expected to lead new insights into 
quantum fundamentals, particularly quantum measurement and quantum-classical transi- 
tions 5]. Recently, Yu and Eberly [6j have showed that two entangled qubits become com- 
pletely disentangled in a finite time under the influence of pure vacuum noise. Zubairy et 
al. [7] have demonstrated how the high quality cavities can be used to realize the new class 
of quantum erasers referred to as quantum disentanglement erasers. Dodd [5] has studied 
the competing effects of environmental noise and interparticle coupling on disentanglement 
by solving the dynamics of two harmonically coupled oscillators. Cucchietti et al. [$] have 
considered the decoherence effect of a non-interacting spin chain on a single qubit. 

In this paper, we study the disentanglement dynamics of a two-qubit quantum sysmtem. 
Here, the key point is that we choose a special correlated XY spin chain to model the 
surrounding environment. This choice of the correlated environment is directly motivated by 
the recent recognition that the single-qubit decoherence induced by a spin-chain environment 



gle-quDit aecom 
ties [9, [lO, 111, [l? 



displays highly interesting properties [9j, llO, llll, due to the unique occurrence of quantum 
phase transition (QPT) in the spin-chain environmental subsystem. Quan et al. [9j] have 
studied the transition dynamics of a quantum two-level system from a pure state to a mixed 
one induced by QPT of the surrounding many-body system. They have shown that the 
decaying behavior of the Loschmidt echo (LE) is best enhanced by QPT of the surrounding 
system. Cucchietti et al. 10] have found that the QPT of the spin-chain environment will 
drive the decay of the quantum coherences in the central quantum system to be Gaussian 
with a width independent of the system-environment coupling strength. 

Motivated by the above-mentioned advances in the QPT effect on the single-qubit deco- 
herence, we turn to study the QPT effect of the environmental spin chain on the two-qubit 
disentanglement of the central quantum system. The coupled spin system we consider in this 
paper consists of two quantum subsystems. One subsystem is characterized by two spin-1/2 
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Hamiltonians, which denotes the general two qubits. We call this subsystem the central sys- 
tem, in the sense that these two spins play the role of measuring disentanglement. Whereas 
the other subsystem (a general XY spin chain in a transverse magnetic field) plays the role 
of the many-body environment. Compared to the Ising model which has been recently used 



to study the QPT effect on the disentanglement [13j, the XY model is parameterized by 7 
and A (see Eq. ([1]) below). Two distinct critical regions appear in parameter space: the 
segment (7, A) = (0, (0, 1)) for the XX spin chain and the critical line A c = 1 for the whole 



family of the XY model 14]. 



The total Hamiltonian for two central spins transversely coupled to a environmental spin 
chain, which is described by the one- dimensional XY model, is given by (h is taken to be 
unity) 

H = - jt (^of ^ + + Aaf) (1) 

N 

1 

= H { E X) + Hj. 

Where given by first line in Eq. ([I]) denotes the Hamiltonian of the environmental 
spin chain, and Hi given by the second line denotes describes the interaction between the 
central two-qubit spins and the spin chain. The Pauli matrices & Z A ^ (a=x, y, z) and of are 
used to describe the central two-qubit spins and the environmental spin-chain subsystems, 
respectively. The parameters A characterizes the intensity of the transverse magnetic field, 
and 7 measures the anisotropy in the in-plane interaction. It is well known that the XY 
spin model described by the first line in Eq. (pQ) encompasses two other well-known spin 
models: the Ising spin chain with 7=1 and the XX chain with 7=0. 
The eigenstates of the operator [o\ + crjj) are simply given by 

\1) = \++)ab, \2) = \--) ab , (2) 
\3) = -j=(\ + -)ab+\ + -)ab), 

\4) = ^={\ + -)ab-\ + -)ab), 

where | ± ±)ab = |±)a ® |±)b denote the eigenstates of the product Pauli spin operator 
°A ® a B w hh eigenvalues ±1. The two-qubit states |1), |2), and |3) are simply spin triplet 



states with total central spin gab = 2, while |4) is singlet state with total central spin 
oab = 0. In terms of these two-spin states, the Hamiltonian ([1]) is rewritten as 

4 

H = J2\j)(j\®Hp\ (3) 

3=1 

where the parameters Xj are 

Ai (2 ) = A ± g, A 3 = A 4 = A, (4) 

and is given from Hjp by the replacement of A with Xj. 

As for quantum criticality in the XY model, there are two universality classes depending 
on the anisotropy 7. The critical features are characterized by a critical exponent v defined 
by £ ~ I A — A c I ^ with £ representing the correlation length. For any value of 7, quantum 
criticality occurs at a critical magnetic field A c =l. For the interval < 7 < 1 the model 
belongs to the Ising universality class characterized by the critical exponent v=l, while for 
7=0 the model belongs to the XX universality class with v—1/2 [l^ . 

Considering the initial state |\l/(O)) = |0s(O)) <8> \iPe(0)), where |0s(O)) is the initial state 
for the two central spins and |^(0)) is the initial state for the environmental spin chain, 
then the subsequent time evolution of the coupled spin system is determined by the time 
evolution operator U(t)=exp(—iHt), \^/(t))=U(t)\^/(0)). Given the central quantity 

for our investigation, i.e., the evolved reduced density matrix for the two central spins, will 
be straightforward to obtain. Thus the key task is to determine the time evolution operator 

n 

in a maximally compact form. For this purpose, we follow the standard procedure [14| by 
defining the conventional Jordan- Wigner (JW) transformation 



of = - 2a+a m ) (ai + a+) , 

1 

Y[(l-2a+a m ) (ai-af), (5) 



m<l 



o\ = -I 



m<l 

h 



of = 1 — 2afai, 



which maps spins to one- dimensional spinless fermions with creation (annihilation) opera- 
tors a~l (o/). After a straightforward derivation, the projected environmental Hamiltonian 
becomes 

N 

H E 3) = -^[( a t+i a i + a t a i+i) + + afat+i) + _ 2a+a;)]. (6) 



Next we introduce Fourier transforms of the fermionic operators described by 
4=t^ J2i a i e ~ i27Tlk/N with k=-M, ..,M and M=(N - l)/2. The Hamiltonian p can be 
diagonalized by transforming the fermion operators to momentum space and then using the 
Bogoliubov transformation. The final result is 



where the energy spectrum £l k is given by 



^ = 2\/(4 X ") 2 + ysin^ (8) 
with ei Aj ^=Aj — cos^, and the corresponding Bogoliubov-transformed fermion operators 



are defined by 

6> (Aj ' 

cos -A^d k — i sin — — d^ k (9) 



with angles satisfying cos9 k —2e k /ft k . It is straightforward to see that the normal 
mode bk,\j dressed by the system-environment interaction is related to the purely environ- 
mental normal mode bk,\ by the following identity 

& M . = (cos 0^)6*,* - i(sma k Xj) )bt k>x , (10) 

where - 9^)/2. 

The time evolution operator for the Hamiltonian (J3J) is then given by 

4 

u® = Y,m\®u%'\t), (ii) 

3=1 

where U^ ,3 \t)=eyjp(—iH^t) is the projected time evolution operator for the spin chain 
dressed by the system-environment interaction parameter Xj. 

Suppose that initially the central spins A and B are entangled with each other but not 
with the spin chain, i.e., at t=0 the two central spins and the environmental spin chain are 
assumed to be described by the product state 

l*t«*(o)) = |0)AB®hM, (12) 



where \4>)ab is the entangled initial state of the two central spins and \^)e) is the initial state 
of the environmental spin chain. The evolved reduced density matrix of the central spins is 
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derived to be 



p AB (t) =Tr E \% ot (t)){% ot (t)\ (13) 

4 

= £ c,4(^i^ (v) (t)4 Aj) wi^)b)0"i, 

where Cj = (7'|0)ab- Equation ({TBI is our starting point for the following derivation and 
discussions. It reveals in Eq. ( Fl3l) that the environmental spin chain only modulates the 
off-diagonal terms of pab through the "decoherence factor" 

F{t) = (i, E \U + E ^'\t)Up\t)\^ E ). (14) 

Whereas, the diagonal terms of pab are not influenced by the environment since for j=j', 
the decoherence factor remains unity. One can see from Eq. ( fl4l) that the decoherence 
factor reflects the overlap between the two states of the environment obtained by evolving 
the initial state with two Hamiltonians and H^ 3 '\ which are different (for j ^ j') 
by the system-dependent parameters Xj and \j> [see Eq. (jjj)]. Furthermore, we notice that 
similar to the single-qubit case, the present decoherence factor F(t) of the two qubits also in 
some special cases has a form of the Loschmidt echo (or fidelity), which can show universal 



behavior (with exponential decay) when are classically chaotic Hamiltonians 15l. Il6| 

The new physical connotation endowed by the special choice of spin-chain environment is 
QPT, which due to its dynamic h ype rsensitivity to the perturbation induced by a single 
qubit as previously investigated {sIIhJ U, 12], or two qubits to be studied here, will play a 



fundamental role in determining the dynamics of the central spin(s) and the corresponding 
decoherence (disentanglement) behaviors. 

Before proceeding the discussion, we would like to point out that the reduced density 
matrix pab sensitively depends through F(t) on the special choice of the initial central-spin 
state \4>)ab and spin-chain state \ip E ). In particular, if \4>)ab lies in the subspace spanned 
by |3) and |4) [see Eq. (j2J)], then there is no dynamic correlation between central spins and 
spin-chain environment, i.e., F(t)=l in this case. Thus we choose the initial state of the 
central spins to have an entangled form 

|0>AB = a|l) + 6|2> (15) 
= a\ + +)ab + b\ )ab- 
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As a consequence, the time evolution of the two central spins will be confined within this 
two-dimensional subspace consisting of |1) and |2) and pab is reduced to a 2 x 2 matrix. 
On the other side, the choice of the initial spin-chain state also needs to be mentioned. 
In the previous work {^], Q, [ll|, 3] involving decoherence of single qubit in the spin-chain 
environment, the qubit is chosen to initially be its unperturbed ground state \g). Then 
\iPe) is naturally and simply chosen to be the ground state of the constrained spin-chain 
Hamiltonian, H g =(g\H\g) . In the present two-qubit case, however, since the initially chosen 
entangled state \4>)ab is not the eigenstate of the unperturbed qubits, thus one cannot choose 
the initial state of the spin chain in the same way as used in the single-qubit discussions. Here 
we choose the initial state \i/je) of the environment to be the ground state \G)\ of the purely 
spin-chain Hamiltonian This choice of \iPe) is natural since it may be assumed that the 
coupling between the central spin subsystem and the spin-chain subsystem is adiabatically 
applied. 

The ground state \G)\ of is the vacuum of the fermionic modes described by 

m ( e w e (x) \ 

bk,x\G)x—0, and can be written as |G}A=lIfc=i ( cos ~2 _ |0)fc|0)-fc + i sin ] , where 

|0)fe and denote the vacuum and single excitation of the kth mode d k , respectively. Note 

that the ground state is a tensor product of states, each lying in the two-dimensional Hilbert 

space spanned by |0)jt|0)_fc and From the relationship between the Bogoliubov 

modes b k ,\ and b ki \. [equation fflQ]) ]. one can see that the ground state \G)\ of the purely 

spin-chain Hamiltonian can be obtained from the ground state \G)\. of the qubit-dressed 

Hamiltonian H E J by the transformation 

l G >^ = IltS 008 ^ +i ^ a k i) K,xfi-k,x j )\G)x r (16) 

Given the initial state |^tot(O)) = |0)AB ® I^e) of the whole system, then our present task 
is to derive the explicit expression for the decoherence factor F(t). First one notices that 
F(t) in Eq. ffMl) can be written as 

\F(t)\ = \x(G\Up X2 \t)U^\t)\G) x \ (17) 
= \x 2 (G\ Yl (cosaf 2) - isma k X2) b_ kM b kM ^ 

k 

x e < 2) V<^n (cos4 Al) +i^\ Xj bt kM ) \G) Xl \. 
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By using the identity e lHxt b^ x e lHxt =b^ x e lUk *, Eq. (17) is rewritten as 

i F wi = IL >0 1 sin ^ Al) sin ^ A2) cos K l} - a i A2) ) ex p H^* + <n M ( 18 ) 

- cosaf sinaf 2) sin (a]j. Al) - a) ) exp (zO Al i + ift A H) 
+ sinaf cosaf 2) sin fa£ Al) - af 2) ) exp (-ift Al t - zfi A2 t) 
+ cosa[ Al) cosaf 2) cos (af - af 2) ) exp (ifij 1 * - ifi Az t) |. 

Equation (ITS]) will be used in the latter discussions, one variant of its form, which will also 
be used for discussion, is the following 



= IW 1 - sin2 K Ai) ) sin2 K i} *) - sin2 H A2) ) sin2 K 2)t ) ( i9 ) 

+ 2 sin (2«[ Al) ) sin ^aj^) sin (fi[ Al) t) sin (of 2) t) cos (tl^k - Sl^t 



4 Al) ) sin (2c4 A2) ) sin 2 U Al) - ) sin 2 (^t\ sin 2 ffl 



fc>0 

Equation ffl9l) [or Eq. (ITS]) ] is one main result in this paper. It can be simplified under 
some special conditions. For example, if one chooses the initial spin-chain state to be 
\iI>r) = \G)\ v then Eq. (TT3]) and corresponding Eq. (fT§|) will be reduced to a LE form given 
in Ref. fll |. It is straightforward to see that each factor Fp. in Eq. (TT9]) has a norm less 
than unity, thus one may expect F(t) to decrease to zero in the large N limit under some 
reasonable conditions. Now we study in detail the critical behavior of the decoherence factor 
F(t) near the critical point A c =l for finite lattice size N of the spin chain. Following Ref. 
9], let us first make a heuristic analysis of the features of F(t). For a cutoff frequency K c 
we define the partial product for F(t) 

\ F M = H^ F k>\nt)i (2°) 

and the corresponding partial sum S(t)=\n \F c (t)\ = - J2£=l l lni? fc|- For sma11 k and sma11 
g (weak coupling) one has 

n[ x) ^2\\-l\+0(k 2 ), (21) 
f^ Aj) ^2\\j - 1| +0(k 2 ), 



s 



and then 




^f2"f7rkg 



(22) 




sm 



iV|(A 1 -l)(A 2 -l)r 



As a result, one has 



S(t) 



^--E(K c ) 1 2 g 2 (X-ir 2 (X 1 -l) 



(23) 



x (A 2 - l) -2 {(A 2 - l) 2 sin 2 (2 |Ai - 1| t) 
+ (Ai-l) 2 sin 2 (2|A 2 -l|t) 



2 |(Ai — 1) (A 2 - 1)| sin (2 |A X — 1| t) 



x sin (2 |A 2 - l\t) cos (4 At)} 



where E (K c ) = Aix 2 K c (K c + 1) (2K C + 1) / (6iV 2 ). In the derivation of the above equation, 
we have omitted the terms related to the sum of k 4 /N A . Consequently, in the short time t 
one has 



when A -»• A c = 1, where r = 8E (K c ) -f 2 g 2 / (A - l) 2 . 

One can see from Eq. ( 1241) that when N is large enough and A — > A c = 1, then |-F C (£)| 
will decay to zero in a short time. It should be noticed that when increasing N, the cutoff 
frequency K c should also linearly increase to remain the validity of Eq. (|24|) . Otherwise, one 
would derive an unphysical conclusion that in the thermodynamic limit, i.e., the number N of 
the sites approaching infinite while keeping the length of the spin chain fixed, r tends to zero 
and thus the approximate expression |.F c (t)| remains unity without any decay. Therefore, 
in using Eq. (1241) to reveal the close relationship between the decaying behavior of |-F(t)| 
and QPT which occur only in the thermodynamic limit, it is necessary to keep the value of 
K c /N invariant when increasing N to infinity. Such kind of scaling relation will be further 
revealed in the latter discussions in this paper. 

Now we check the dynamical property of |-F(t)| by numerical analysis calculated from 
the exact expression Eq. (f!9l) . In Fig. 1(a), the |-F(t)| is plotted as a function of magnetic 
intensity A and time t for iV=201, g=0.05, and 7=1.0 (i.e., the case of Ising spin chain 



F c (t)\ 



e 



(24) 
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FIG. 1: (Color online), (a) Disentanglement factor |-F(£)| as a function of magnetic intensity A and 
time t for two central spin qubits coupled (with coupling strength g = 0.05) to an Ising (7 = 1.0) 
spin chain with the size N = 201. (b) Disentanglement factor for different sizes of Ising spin chain 
at QPT point (A = 1). 

and in the weak coupling regime). One can see that apart from the critical point A c , the 
|.F(t)| in time domain is characterized by an oscillatory localization behavior. When the 
amplitude of A approaches to A c , the degree of localization of |-F(t)| is decreased to zero. 
The fundamental change occurs at a critical point of QPT, ie., A=A C =1. At this point, as 
revealed in Fig. 1(a), the |-F(t)| evolves from unity to zero in a very short time, which implies 
that the disentanglement of two central spins is best enhanced by QPT in the environmental 
spin chain. The size dependence of the decoherence factor is shown in Fig. 1(b) for A=A C 
and g=0.05. Not surprisingly, with increasing iV towards thermodynamic limit, the role of 
QPT in Ising spin chain becomes clear by completely disentangling the two central qubits 
in a very short time. 

As mentioned at the beginning of this paper, for the XY model we employed, there are 
two distinct critical regions in parameter space. Region I is the segment (7, A) = (0, (0, 1)) 
for the XX spin chain, while region II is a critical line A c =l for the whole family of the XY 
model (including the special case 7=1 of Ising model). We find that the best-enhancement 
behavior of the disentanglement factor |-F(t)| only occurs in QPT region II except for the 
point (0,1). Whereas in the whole region I and at the point (0,1), |-F(t)| remains unity 
during the time evolution, and thus the QPT in the environmental spin chain has no any 
effect on the entanglement of the two central spins. This full localization behavior of |-F(t)| 
can be seen from the analytic expression, Eq. (|2^|) . in which r=0 for 7=0, indicating no 
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FIG. 2: (Color online). Disentanglement factor |-F(t)| as a function of spin anisotropy parameter 
7 and time t for two central spin qubits coupled to an XY spin chain. The other parameters are 
set to be A = 1, N = 201, and g = 0.05. 

decay in regardless of the variation of A and the couping strength g. Physically, this 

vanishing of disentanglement for the two central qubits under an XX spin-chain environment 
can be seen by noticing that the parameters 0^ and 9^ in Eq. fflOl) are zero (or it) at 
7=0. In this case, the fermionic modes bf.^. and b^x coincide each other, which leads 
to complete overlap between the ground state \G)\ j of and the ground state \G)\., of 
H^ J '\ \G)\ j =\G)\. l . As a result, one sees from Eq. (THI) that the disentanglement factor F(t) 
keeps an invariant value of unity during its time evolution. Thus one arrives an important 
conclusion that the enhancement of the disentanglement by QPT may be broken by special 
choice of the spin-chain the occurrence of ground-state "accidental" degeneracy among the 
system-dressed environmental spin-chain Hamiltonians in critical parameter space. 

For further illustration, we show in Fig. 2 \F(t)\ as a function of time and 7 for A=1.0, 
iV=2GT, and g=0.05 (weak coupling), which corresponds to critical region II. One can see 
that with deviating 7 from zero, the disentanglement factor gradually evolves towards zero 
in an oscillatory way. 

After discussing the QPT effect on the disentanglement of two spin qubits in weak cou- 
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400 



FIG. 3: (Color online). Disentanglement factor |-F(i)| as a function of time in strong coupling 
regime. The system parameters are chosen to be A = 1.0, N = 201, 7 = 1.0, and g = 500.0. The 
exact numerical result is shown by solid line, while the approximate Gaussian envelope factor is 
plotted by dashed line. 



pling regime (g << A c ), we turn now to study the QPT effect in strong coupling regime 
(g » A c ). In Fig. 3 (solid line) we display the time evolution of |-F(£)| for the values of 
A=1.0, 7=1.0 (Ising model), iV=201, and g=500. Besides the best-enhancement behavior 
(|.F(t)| — ^0 in final time) of the disentanglement as discussed above, one additional promi- 
nent new feature, which is absent in the weak coupling case, is that the decay of |-^(t)| is 
now characterized by an oscillatory Gaussian envelope. To explain this, we starts from the 
observation that when g 3> 1, the spin-chain energy spectrum in Eq. ([8]) can be simplified 



to n[ Xl) » 2e[ Xl) and fi£ Aa) w -2ef 2 } . Thus from Eq. Q one has 9 ( k Xl) 



and Of 2) « 7T. 



This leads to the approximate identity a 
Eq. (TI5I) one can obtain 



(Ai) 
k 



a 



(Aa) 



-7r/2, by substitution of which into 



F (t)\ « JJ I cos 2 4 Al) exp (zO fc t) + sin 2 a£ l) exp (-ifi*,*) 



(25) 



f^ 1 +^fc 2 . Remarkably, the above expression for \F (t)\ is completely analogous 



)A 2 



where flk 

to the one found when studying deco 
(see Eq. (16) in Ref. 



environment 
derivation given in Ref. 
is as follows 



lerence on a qubit induced by noninteracting spin 
). Thus, one can exactly follow the mathematical 
and Ref. [10j. The resultant approximate expression for 



\F(t) \ « exp (s 2 N t 2 /2) |cos(fi*)| 



(N-l)/2 



(26) 
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where Q is the mean value of Qk, i-e., £l=-h J2k>o^ k ' anc ^ 

4 = E sin22 4 Al) ^ 2 - (27) 

fc>0 

Here the quantity Sk describes the deviation of from its mean value Q. It is straightforward 
to obtain Q m 4g + 7 2 / 'g and 5k ~ — ^-cos^. We remark that the present Gaussian 
character of the disentanglement factor is not only confined to the QPT regime. Here it is 
mainly for purpose of the consistency in organizing this paper that we focus our attention to 
the strongly coupling behavior of |-F (t) | in the vicinity of QPT. After a careful analysis of Eq. 
([27)1 . we further find that the width of the Gaussian envelope is proportional to g , y~ 2 N~ 1 / 2 , 
which is an important scaling relation between the decaying factor |.F(t)| and the system 
parameters in the strong coupling QPT regime. For comparison with the exact numerical 
result, we also show in Fig. 3 (dashed line) the Gaussian envelope factor exp (— s 2 N t 2 /2) in 
the approximate expression (1261) of |F(t)|. Clearly, the agreement is very good, indicating 
the validity of our approximation in QPT region (A=l for Ising model) with strong system- 
environment interaction. Figures 4(a)-(c) display the exactly numerical results (solid lines) 
of and the analytic results (dashed lines) of Gaussian envelope factor exp (—s%t 2 /2) 

for A=l and different choices of the other parameters g, 7, and N. It remarkably reveals in 
Fig. 4 that the decaying width of is proportional to the product g'j^N' 1 ^ 2 , exactly 

as we have analyzed in the above discussions. 

Finally, we find that in the vicinity of QPT, the shape of the disentanglement factor 
|F(t)| during its time evolution is invariant under the scaling transformation t — > t/a, 
5 — > aS, g — > oeg, and 7/iV — > acy/N, where 5 = X c — A characterizes the vicinity of 
QPT. To illustrate this remarkable scaling property, we plot in Figs. 5 the exact numerical 
results of evolution of for different values of the system parameters. Here the values 

of the system parameters used in Fig. 5(b) are obtained from those used in Fig. 5(a) by a 
scaling factor a = 0.1. Clearly, it shows in Fig. 5 that the exact time evolution of |i^(i)| 
faithfully follow this scaling transformation. Remarkably, the similar scaling property has 
been recently found [9] in studying dynamics of the LE for a single qubit coupled to an 
Ising-type spin chain. Clearly, this scaling rule in the disentanglement factor for two 

entangled qubits or in the LE for the single qubit is highly meaningful in quantum computing 
and quantum information processing. 

To understand this scaling property, here we give a detailed analysis of the behavior of 
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FIG. 4: (Color online). Disentanglement factor |-F(i)| as a function of time in strong coupling QPT 
regime (A = 1.0) for different choices of parameters 7, g and N to show their relationships with 
the decaying width of |-F(i)|. Again, the exact numerical result is shown by solid line, while the 
approximate Gaussian envelope factor is plotted by dashed line. 

|.F(t)| in the vicinity of the critical point A c =l in the case of weak coupling strength g. Note 
that although in the present context we only concern the specific model employed in this 
paper, the following analysis can be easily applied to the other cases. We first notice that 
in the expression of |-F(t)| [Eq. ( fl9l) ]. most factors remains nearly unity. Thus only very 
few FkS have remarkable effect on the shape and amplitude of |-F(t)|. From Eq ( fl9l) one 
can see that in order for the factor to deviate prominently from unity, at least one of its 
two coefficients sin2ajj. Aj ' ) (j = 1,2) should be considerably non-zero. Next let us check the 
value of sm2a^ 3 \ For this we define k^ 3 ^ which enables = Aj — cos \ 2irki X ^ /n\ as 
small as possible. For small 5 (i.e., A — A c ) and g, one can see that ki X ^ <C M. From the 
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FIG. 5: (Color online). Scaling behavior of 1-^(^)1 in the vicinity of the critical point A c =l in 
the weak coupling case. The parameters used in plotting curves in (b) are related to those used 
in plotting curves (with the same curve type) in (a) by the transformation g — > ag, S — > a.5, 
j/N — > aj/N with a = 0.1. One can see that by further transformation t — > t/a, figures (a) and 
(b) will completely overlap. 

definition of , we can write down 

sin 2a[ Aj) = T^g sin {2-nk/N) (28) 

for j — 1, 2, respectively. One can see from Eq. (128]) and the expressions of Q { k Xj) and 
Q k that for small g, to enable sm2a^ considerably non-zero, three conditions should be 
satisfied: (i) k should be close to kc and k^ in order for the amplitude of 7 sin (2nk/N) 
to be comparable with ej^ and ej^; (ii) g is small enough so that k could be close to ki X ^ 
and kc at the same time, (iii) 5 is small which leads to small value of sin (2irk/N) when k 
approaching ki and ki X \ Under these three conditions, one has the following approximate 
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expressions 

nJM « 2 [(5 T g) 2 + 4 7 VA; 2 /iV 2 ] 1/2 , (29) 
fi^2(5 2 + 4 7 2 7r 2 A; 2 /iV 2 ) 1/2 . 

Combining Eq. (1251) and Eq. (1291) . one immediately finds that the transformation g — > ctg, 
5 -> a<J, and 7 /iV -> a 7 /iV leads to fi£ Aj) -> orfij^, fi£ -> aft£, while sin2o4 Aj) and 
cos2a^ 3 ^ remaining invariant. As a result, the time evolution of in Eq. (11 9p is well 

invariant under further transformation t —>■ t/a. This is what one has seen from the exact 
results in Fig. 5. 

In summary, we have studied the dynamic process of the disentanglement of a coupled 
system consisting of two spin qubits and a general XY spin chain. The exact expression 
of the disentanglement factor |.F(i)| has been obtained. The relation between |.F(t)| and 
the QPT in the environmental spin chain has been extensively illustrated. It has been 
shown that in general, the disentanglement of the two qubits is best enhanced when the 
environmental spin chain is exposed to QPT in either strong or weak coupling case. Both 
the heuristic analysis and numerical calculations have shown the sharply decaying behavior 
of the decoherence factor in the vicinity of the critical line A=A C =1. This decaying behavior, 
on the other side, has been found to break for the particular XX spin chain ( 7 = 0), in which 
case is not influenced by the environment. In the strong coupling case, it has been 

numerically and analytically found that in the vicinity of QPT the disentanglement factor 
decays to zero in an oscillatory Gaussian envelope. The width of the Gaussian envelope has 
been found to scale with a form g^~ 2 N~ l l 2 . Furthermore, we have established a scaling rule 
for the time evolution of the disentanglement factor in the vicinity of QPT. We expect that 
the present results may shed light on the role of strongly correlated environment played in 
the disentanglement dynamics of multi-qubits. 
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